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Abstract 

We present basic properties of Dipolar SLEs, a new version of 
stochastic Loewner evolutions (SLE) in which the critical interfaces 
end randomly on an interval of the boundary of a planar domain. We 
present a general argument explaining why correlation functions of 
models of statistical mechanics are expected to be martingales and we 
give a relation between dipolar SLEs and CFTs. We compute SLE 
excursion and/or visiting probabilities, including the probability for 
a point to be on the left/right of the SLE trace or that to be inside 
the SLE hull. These functions, which turn out to be harmonic, have a 
simple CFT interpretation. We also present numerical simulations of 
the ferromagnetic Ising interface that confirm both the probabilistic 
approach and the CFT mapping. 



Two types of stochastic Loewner evolutions (SLE) have been thoroughly 
studied, the chordal and the radial SLEs [H El El Ej . The former describe ran- 
dom curves joining two points on the boundary of a simply connected planar 
domain while the latter describe curves joining one point on the boundary 
to a point in the bulk of the domain. They correspond to two inequivalent 
normalizations of conformal maps between simply connected domains in C. 
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Using geometrical constraints, we realized in ref . jH] that there was a third in- 
equivalent normalization and thus that there was yet another process, which 
we called the dipolar SLE, with all properties required to define an SLE, 
depending as usual on a real positive parameter k. Concretely, dipolar SLEs 
describe random curves in a simply connected planar domain which start 
at a point x on the boundary and are stopped the first time they hit the 
boundary on a specified interval not containing xq. They generalize 

chordal SLEs. One of the virtues of dipolar versus chordal SLEs is that the 
hull does not fill the full domain for k > 4. This makes the dipolar geome- 
try physically appealing. It is for instance the most natural one to describe 
several interface properties already computed, such as Cardy's formula for 
percolation jHUH]- The chordal case corresponds to the limit when x + and x_ 
merge together, and in this limit the hull invades the full domain for k > 4. 

The aim of the following is to introduce dipolar SLEs and to present their 
basic properties. 

Section [T] gives the precise definition of dipolar SLEs. 

Section |21 is a digression of general nature which emphasizes why (con- 
ditional) correlation functions of models of statistical mechanics are martin- 
gales for appropriate stochastic processes. This can be used to elucidate the 
link between SLEs and conformal field theories (CFTs). 

Section El applies these general ideas to build the CFT interpretation of 
dipolar SLEs, identifying in particular the boundary operators acting at x 
and x±. 

Section |U is devoted to the computation of some basic dipolar bulk prob- 
abilities for k > 4, namely the probability for a point to be on the left of 
(resp. on the right of) or swallowed by the SLE hull. These probabilities 
are computable because they turn out to be harmonic solutions of the gen- 
eral martingale equation (|7J). Why this equation has interesting harmonic 
solutions remains to be explained. In fact, the probabilities for a point to 
be swallowed by the hull from the left or from the right for k > 4 are non- 
harmonic, and so are the probabilities to be on the left or on the right of 
the SLE hull for k < 4. The explicit computation of these probabilities has 
eluded us. 

Section is devoted to the limiting case k = 4 and its relation to free 
field theory with an alternation of appropriate boundary conditions, namely 
Dirichlet between X- and xo, and between x and x + (but with two different 
boundary values), and Neumann between x + and x_. 

In Section |U] we compute, for arbitrary k, the distribution of the hitting 
point of the hull on the interval that does not contain xq. We also 

give a detailed CFT derivation of the result. 

In Section Owe present Monte Carlo simulations of the Ising ferromagnet 
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with specific boundary conditions corresponding to the case k = 3. The 
distribution of the interface endpoint agrees very well with the theoretical 
prediction presented in this article. This confirms the validity of the mapping 
to the Ising model. 



1 Dipolar SLEs. 

As for any SLEs, the random curves, often called SLE traces, are encoded into 
a family of conformal maps gt(z) parameterized by a 'time' t. When the traces 
are simple curves, the map gt uniformizes the complement of the portion 
of the curve 7(0^1 in the domain on which the curve is growing back into 
this domain. Such maps exist by the Riemann mapping theorem and they 
are defined up to SL 2 (M.) transformations by global conformal symmetry. 
Thus we need three conditions to fully specify them. For dipolar SLEs, we 
choose two points on the boundary, x±, and impose that gt(x±) = x± and 
g' t (x+) = g' t (x-) for each t. The process is then defined by specifying the 
stochastic evolution of g t via a Loewner equation. The dipolar SLE maps 
describe curves starting at a point x 7^ x± on the boundary and ending on 
the boundary at a random point on the interval between the two fixed points 
that does not contain Xq. The crucial point is that for a simply connected 
domain, the group of conformal automorphisms fixing two boundary points is 
a one parameter Lie group isomorphic to the additive group of real numbers. 
Hence, there is a canonical definition of a Brownian motion between x_ and 
x + starting at Xo, at least when the corresponding boundary interval is a 
simple curve. See jH El El Hj for basic - and not so basic - material on SLEs. 

To be more concrete, let §a = {z G C, < Smz < 7rA} be the strip of 
width 7rA which is the geometry adapted to dipolar SLEs. By definition, the 
dipolar Loewner equation in §a reads: 

d t gt(z) = - — cWoAV gt=o{z) = z, (1) 

tanh((0 t (z) -&J/2A) 

with £ t = y/HB t with B t a normalized Brownian motion and k a real positive 
parameter so that E[£t£ a ] = Kmin(t, s). The two boundary fixed points are 
x± = ±00 and the starting point xq is the origin. The maps gt are normalized 
to fix x±. For fixed z, gt{z) is well-defined up to the time r z for which 
g Tz (z) = C, Tz . Times r z are called swallowing times. In the limit A — > 00 we 
recover the chordal SLEs. A is simply a dilatation factor. Unless otherwise 
specified, we will set A = 1 in the following and we will look at dipolar SLEs 
in the strip Si. 
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As for chordal or radial SLEs [U 121 EH E] , the SLE hull is reconstructed 
from g t by K t = {z such that r z < t} and the SLE trace 7[ ,t] by 7(2) = 
lim e ^ + (£ t + ie). It is known that 7[o,t] is almost surely a curve. It is non- 
self intersecting and it coincides with Kt for < k < 4, while for 4 < k < 8 
it possesses double-points and it does not coincide with K t which is then the 
set of points swallowed up to time t. For k > 8, the trace is space filling. 
However, in contrast to the chordal SLE case, for dipolar SLEs, the hull does 
not fill the entire domain even for k > 8. 

It will be convenient to consider the map which maps the tip of the SLE 
trace back to its starting point. Thus we translate g t and let ft(z) = gt{z)—^ t 
so that ft(j(t)) = 0. The dipolar stochastic equation is simply: 

df t (z) = A' 1 coth(/ t (z)/2A) + d&. 

The maps f t are such that, for s > t, f s o ff 1 is distributed as / s _ t . 

Eq.([T} may be integrated in the simple deterministic case with £t constant. 
Then cosh(g t (z) — £)/2 = e'/ 2 cosh(z — £)/2 for A = 1. The trace j(t) is 
determined by cosh(7(t) — ^)/2 = e~'/ 2 , so that it starts on the real axis at 
£ and stops when it touches the upper boundary of the strip at the point 
in + £. As expected, all points at the left of the trace are mapped to the 
fixed point at — oo by f t as t — > oo, while those on the right of the curve are 
mapped into the other fixed point +oo. 

Let us return to the case when £ t is a Brownian motion. A simple proba- 
bilistic argument, explained in [Hj, shows that the probabilities for the trace 
to touch the upper boundary in finite time vanishes. So, the trace only 
touches this boundary at infinite time. This amounts to stop the process 
when the trace hits the upper boundary, a criterion which is invariant under 
time reparametrizations of the evolution. Since locally the dipolar SLE trace 
looks the same as the chordal or radial SLE traces this leads to the fol- 
lowing picture. For k < 4, the hull coincides with the trace. The latter does 
not touch the lower boundary, ie. the real axis R, but stops when it hits the 
upper boundary, ie. in + R, at same random point. For k > 4, the trace 
is not a simple curve and points in Si are swallowed, the set of which form 
the SLE hull. The hull then intersects the lower boundary, ie. the real axis, 
and the trace hits this boundary an infinite number of times but hits only 
once the upper boundary, again at some random point, and stops there. See 
Figure (JTJ). 
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Figure 1: Two schematic examples of dipolar SLEs. 



2 Statistical mechanics and martingales. 

This section is a recreative interlude in which we explain why (conditional) 
correlation functions of models of statistical mechanics are martingales for 
appropriate stochastic processes. As a very general statement this remark 
may look tautological but it is nevertheless quite instructive. In particular 
it provides a key to decipher the relation between SLEs and CFTs. 

Let C be the configuration space of a statistical model. For simplicity we 
assume C to be discrete and finite but as large as desired. Let W c be the 
Boltzmann weights and Z the partition function, Z = Y2cec W c - 

We imagine having introduced a family Q t of partitions of the configura- 
tion space whose elements C at are labeled by indices a t : 



The index t, which will be identified with 'time', labels the partitions. By 
convention Qq is the trivial partition with C as its single piece. We assume 
these partitions to be finer as t increases, which means that for any s > t 
and any element C at of the partition at time t there exist elements of Q s 
which form a partition of C at . An example of such partitions in case of spin 
statistical models consists in specifying the values of local spin variables at 
an increasing number of lattice points. Block spin clustering used in renor- 
malization group is another way to produce such partitions. A SLE inspired 
example consists in specifying the shapes and the positions of interfaces of 
increasing lengths. 



C 




C at disjoints. 
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We define the restricted partition function Z at by 

Since restricting the summation to a subset amounts to impose some condi- 
tion on the statistical configurations, Z at is the partition function conditioned 
by the knowledge specified by C at ■ 

To define a stochastic process we have to specify the probability space and 
a filtration on it. The simplest choice is C as a probability space equipped 
with its canonical er-algebra, ie. the one generated by all its singletons, 
and with the probability measure induced by the Boltzmann weights, ie. 
P[{c}] = W c /Z. In particular the probability of the event C at is the ratio of 
the partition functions 

P|A*t] = z a jz. 

To any partition Q t is associated a a-algebra on C, ie. the one generated by 
the elements of this partition. Since these partitions are finer as 'time' t in- 
creases, it induces a filtration Tt on C equipped with its probability measure. 

Now, given an observable O of the statistical model, ie. a function c — > O c 
on the configuration space, we can define its conditional average 

(0) t = E[0\T t }. 

By construction, (0) t is a function which is constant on any element C at of 
the partition Q t and takes values 

(0)t\ Cat = -r- E °° W - 

This is simply the statistical average conditioned on the knowledge specified 
by C at . 

By construction, (0) t is a (closed) martingale with respect to T t - Indeed, 
for t > s, 

mo) t \F s ] = E\E[0\F t ]\r 8 ] = E[0\F 3 ] = (0) s , 

where we used standard properties of conditional expectations and that 
T t C T s for t > s. In particular, one may verify that its average is time 
independent and equals to the statistical average: 

E[(0> t ] = EPfCcJ^k, 

at 

= ^O c W e = (0). 

cec 



6 



This observation formally applies to critical interfaces and hence to SLEs. 
The remarkable observation made by O. Schramm is that conformal invari- 
ance implies that the filtration associated to the partial knowledge of the 
interface is that of a continuous martingale, i.e. that of a Brownian motion 
if time is chosen cleverly. The only parameter is k. The physical parameters 
of the CFT, for instance the central charge, can be retrieved by imposing 
that the correlation functions (0) t conditioned by the knowledge of J- t be 
martingales. The CFT situation is particularly favorable in that going from 
(0)t to (O) is pure kinematics. 



3 CFT connections. 

Connection with conformal field theories can be done using the CFT operator 
formalism. The latter is simpler if the boundary CFT is considered in the 
upper half plane H = {z e C, Qmz > 0}. 

SLE processes are transported onto any simply connected domain by 
conformal transformations, by definition. In H, we fix the boundary points 
x± to be ±1 and Xq to be the origin. The conformal map uniformizing 
§a onto EI is <f(w) = tanh(w/2A). Let gt be the dipolar SLE map in H, 
gt o ip = ip o g t . The stochastic Loewner equation for dipolar SLE in HI reads 

f l-g t (z) 2 \ f l-g t {z)t*nht t /2 \ „ 

The maps gt are normalized to fix x± = ±1 and to have equal derivatives at 
these two points : $(±1) = e~*. 

The map ft can also be transported in the upper half plane to produce a 
map ft = tanh(/ t /2A), which also fixes x± = ±1 but which maps the tip of 
the SLE trace in H to the origin. Its expression is: 

if s 9t(z) -tanh^/2 
Jt[ z ) 



l-gt(z) tanh^/2" 



The stochastic differential equation that f t satisfies directly follows from that 
obeyed by f t . 

As for the chordal and radial cases [Zj, the connection between SLEs and 
CFTs may be established by associating to ft an operator Gt which imple- 
ments this conformal transformation in the CFT Hilbert space. Since f t fixes 
the point x + , or x„, G t may be constructed as an element of the enveloping 
algebra of the appropriate Borel sub-algebra of the Virasoro algebra. See [Z| 
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for details. It intertwines between primary fields, say z) of conformal 

dimensions (h, h), and their images under f t . Namely, 



G^<S> h>n (z,z) G t = [fl(z)] h [fl(zW **,&(/*(*),/*(*))■ 

The stochastic differential equation that Gt satisfies directly follows from 
that of ft, or g t . It reads [Bj: 

G" 1 dG t = (-2W- 2 + ^W\)dt - W-idtt, (2) 
where W-2 and W-\ are elements of the Virasoro algebra 5 , 

= \(L- 2 - Lq), W- X = i(L_! - L x ). 

In this algebraic setting, the differences between radial and dipolar SLEs may 
be viewed as coming from a different choice of real forms in the Virasoro 
algebra. 

As in the chordal and radial SLEs, the key point is now the construc- 
tion of a generating function of local martingales which is obtained using a 
representation of the Virasoro algebra degenerate at level two. 

Let \u) be the highest weight vector of the irreducible Virasoro module 
with central charge c and conformal weight hip, 

(k-6)(8-3«) , 6-k 

C = , Ail -o = , 

2k ' ' 2k ' 

then 

M t = e +2h °-^ t G t \u) (3) 
is a local martingale, with 

(6 -«)(«- 2) 



ho-,1/2 



16k 



This follows from the null vector relation (— 2L-2+^L^)\u) = 0. Indeed, 
a simple rearrangement leads to 



2 -1 = i(" 2L -2 + 2 L ' i) + + 8 (L ' " 2L - lLl ^' 



We denote by L n the generators of the Virasoro algebra with commutation relations 
— ? 

12 1 



[L n , L m ] — (n — m)L rn+n + -^n(n 2 — l)<5„+ m ,o and c is the central charge 
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so that (-2W^ 2 + § W 2 _ 



i 



2/i 0;1 / 2 |a;) and 



G -1 dG t |w) = -d*2/io;i/2|w) + dCtW_i|w). 



In particular, by projecting this local martingale on vectors (v\ and as- 
suming appropriate boundedness conditions, we get that the expectations 
E[ e 2fc. ;i/2* (vjG'^a;}] are time independent. 

The prefactor e +2fe ° ;1 / 2 * accounts for the insertion of two boundary con- 
formal fields, each of dimension foo;i/2, localized at the two fixed points x± 
since f' t {x + )f' t {x-) = e~ 2t . Alternatively, the local martingales can be writ- 
ten as the ratio of two correlation functions with insertions of the boundary 
operator ^ 1;2 creating the state \cu) at x and of the operators ^051/2 at the 
fixed points x±: 



for any operator O. Since martingales are key ingredients for computing 
probabilities, this statement implies that dipolar SLE probabilities will be 
expressible in terms of CFT correlation functions with insertions of these 
boundary operators. 

By conformal invariance, this last expression can be transported to any 
simply connected planar domain §. If g t is the conformal map describing the 
growth of the dipolar SLE hulls in §, mapping the tip of the SLE trace to 
xq G <9S>, then the ratios 



are local martingales. Here ^0*1/2' ^1*2 or ® 9t denote the images of the corre- 
sponding fields by g t , ie. their pull-back by g t . As explained in the previous 
section, this statement has a natural explanation in basic statistical mechan- 
ics. Although this construction can be seen merely as a trick to construct 
martingales, the previous section explains why the martingale property is 
what ensures that the CFT identified in this way is precisely the critical con- 
tinuum limit of the statistical mechanics model that produces the interface 
described by SLE. 

— Examples: 

Dipolar SLEs has a simple statistical interpretation for k = 3 which cor- 
responds to the Ising model with c = 1/2. Then ipi-2 with hi- 2 = 1/2 is 
the boundary condition changing operator between spin (+) and spin (— ) 



(X_ )V>Q;i/2 (x+) O G t ipxp (x )) 
{^0;l/2(X-)lp .i/ 2 (x + ) G t Iplflixo)) 



(4) 



{^i /2 (x-)^i /2 (x + )Q 9 ^iA^)) 

(^0-1/2 ( X ~ )^0;l/2 ( X + )^1;2 (x ) ) 



(5) 



9 



fixed boundary conditions, while i/)q-i/2 with ^o ; i/2 = 1/16 is the boundary 
condition changing operator between fixed and free boundary conditions [9 . 
Hence, along the boundary of the domain one encounters the boundary con- 
ditions fixed (+), then free, and then again fixed but (— ), as depicted in 
Fig-O- It is thus clear that we need one operator ^i ;2 and two operators 
V ; 0;i/2 to describe this system. 

As we shall see in the following sections, the = 4 which corresponds 

to a free bosonic field with c = 1 also has a simple interpretation with an 
alternation of Dirichlet and Neumann boundary conditions. 

For the 3-state Potts model we may propose the following interpretation. 
The microscopic spin variables take three possible values (1), (2) and (3) 
related by 5*3 symmetry. The 3-state Potts model corresponds to k — 10/3 
with central charge c = 4/5. Then /i 1;2 = 2/5 and /io ; i/2 = 1/15. The 
operator ipi- 2 is the boundary operator between a fixed boundary condition 
with all spins (1) and a mixed boundary condition (2 + 3) with a mixture of 
spins (2) and (3). The operator ^o;i/2 is the boundary operator between the 
boundary conditions (1) and (1 + 2). It is also the lowest boundary primary 
operator generated by a change of mixed boundary conditions from (1 + 2) 
to (2 + 3), see ref. 9j. Thus the dipolar SLE with k = 10/3 should correspond 
to the 3- Potts models with the succession of boundary conditions, fixed (1), 
mixed (2 + 3) and mixed (1 + 2). 

Unitary minimal CFTs with c = 1— 6/m(m+l) with m integer correspond 
to two values of k related by duality: k = 4(m + l)/m or k, = Am/ (m + 1). 
Since the identification of ^o ; i/2 depend on the parity of m — ie. /io ; i/2 = 
h m /2-m/2 for m even, k > 4, but h . 1/2 = h( m +i)/2;(m+i)/2 for rn odd, k < 4, - 
a simple microscopic interpretation is possible only for one of the two choices 
of K. For instance, although k = 16/3 corresponds to the Ising model, the 
role of the Virasoro representation with weight /to;i/2 — 5/192 in the Ising 
model has not yet been clearly identified. 

4 Bulk visiting probabilities and harmonic func- 
tions (k > 4). 

Let us now look at SLE bulk properties. We assume n > 4 and we deal with 
dipolar SLEs in the strip Si. We shall evaluate the following probabilities: 

(i) The probabilities Pi(z,z) (resp. P r (z,z)) for a bulk point z not to 
be swallowed by the SLE trace and to be on the left (resp. the right) of 
the trace. This is the probability for the point z to be on the left (resp. 
the right) of the exterior frontier of the SLE hull viewed from the boundary 
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point a;_ (resp. x + ). It is therefore the probability for the existence of a path 
joining x_ (resp. x + ) to the boundary interval leaving the point 

z on its right (resp. left) and included into one cluster of the underlying 
model of statistical mechanics. For k = 6 it bears some similarities with 
the probability computed by Smirnov |Hj to prove the equivalence between 
critical percolation and SLE & . 

(ii) The probability Pi n (z, z) for the point z to be in the SLE hull. We do 
not distinguish the events in which the point has been swallowed from the 
right or from the left. Since this probability is also that for the point z to be 
in left or right frontiers of the hull it gives informations on the shape of the 
hull. These probabilities turn out to be harmonic functions for all values of 
k > 4 and are proportional to the CFT correlation functions 

(tpO;l/2(x_)llJo;l/2(x + )^o(z, z)lfafl(x )) (6) 

involving a weight zero bulk primary field $o- 

As usual, a way to compute these probabilities is to notice that the process 
t — > P(f t (z), ft(z)) is a local martingale. Indeed, since f s o ff , s > t, is 
independent of f t and distributed as f s - t , the function P(f t (z), ft{z)) is the 
probability of the events (i), or (ii), conditioned on the process up to time t 
and, as such, it is a martingale. As a consequence, the drift term in the Ito 
derivative of P(f t (z), f t (z)) vanishes which implies that P(z,z) satisfies the 
following differential equation: 

Kd,8 g P + (coth | + ^d z )d z p + (coth z - + ^d 2 )d- z P = 0. (7) 

The main observation in this section is that, quite remarkably, eq.((ZJ) has 
interesting harmonic solutions, in fact enough harmonic solutions to compute 
Pi, P r and P in . 

The boundary conditions depend on which probabilities we are comput- 
ing: 

(i) For the probability to be on the left of the hull, this requires: 

P,(-oo) = l, P(+oo)=0, F,(0) = 0. (8) 

Similar conditions hold for P r (z, z). 

(ii) For the probability to be in the hull, it requires: 

P in {-oo) = P in (+oo) = 0, P in (0) = 1. (9) 

These boundary conditions follow by noticing that if point z is swallowed 
at time r z then \\m. t / Tz ft(z) = 0, if it is not swallowed but is on the left 
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of the trace then lim^oo ft{z) = — oo, and if it is not swallowed but is on 
the right of the trace then hm t y 0O / t (2;) = +00. These conditions are such 
that at the stopping t z = min(r z , 00) the martingale P(ft(z), ft(z)) projects 
on the appropriate events (i) or (ii), ie. P(f? z (z), fr z (z)) = lcvents- As a 
consequence, the probability of these events are: 

P[events] = E[l cvcnts ] = E[P(f ?z (z)J % (z))\ = P(z,z). 

The martingale property has been used in the last equality 

— Case (i): probability to be on the left of the trace. 
The solution of the martingale equation (|7J) satisfying the boundary con- 
ditions (jHJ) is given by the harmonic function 

p,M) = i-|^ (io) 

xsmh (00) 

with 

\ f z du 
F z) : 



^ (sinhu/2) 4 / K ' 

The function F(z) is well-defined and analytic on the strip §! for all k's. For 
n > 4, F(z) is bounded and has a continuous extension to the closure of 
Si. For k < 4, it is unbounded near the origin and so are the corresponding 
solutions to eq.((ZJ). In that case these solutions only lead to local martingales 
and not true martingales. See Appendix for details on its definition and on 
its properties. We have QmF(oo) = — sin(27r//t)) I = — sin(47r/K)) J with 

/ + OO /"OO 
dy(coshy/2)- 4 / K , J= dy(smhy/2)-^ K . 
■00 Jo 

As a check one may verify that Pi(z, z) behaves as expected on the boundary. 
On the positive real axis, (sinh z/2) is real and positive so that 

Pi{x) =0, x e R+, 

in accordance with the fact that no point on the real axis can be on the left 
of the trace. On the negative real axis, (sinhx/2) 4 / K = e i47r / K (sinh |x|/2) 4//k 
and 

1 f +co dy 

which interpolates between one and zero. It gives the probability of the 
hull not to spread further than x on the negative real axis. On the upper 
boundary, 

PAz = itt + x) = 1 - - / , * , zetn + R, 



I (coshy/2) 4 /- 
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since there (sinh z/2) 4 / K = e 2Mr//K (coshx/2) 4 / K . This yields the density proba- 
bility for the trace to stops on an interval [x, x + dx] on the upper boundary. 

— Case (ii): probability to be inside the hull. 

The solution of the martingale equation (J7J) satisfying the boundary con- 
ditions is given by the harmonic function 

P (z z) - ^^IIMJ 
in[ ' ] ~ 3ta[e«*/"F(0)] 1 ] 

with identical function F(z) as above and ^sm[e l2n ^ K P(0) ] = — sm(2n/K) J. 
Again, Pj n has the expected behavior on the boundary. Since e t27T ^ K F(z) is 
real on the upper boundary, we have 

P in (z, z) = 0, z E in + R, 

in agreement with the fact that no point on the upper boundary can be 
swallowed. Pj n is even on the real axis and 

1 r +oc dv 

This is of course complementary to Pi(x) for x negative. 

— CFT interpretation. 

The correlation function (JHJ) has a natural interpretation in the Coulomb 
gas representation of CFTs: the weight zero primary field $o is simply the 
integral of the screening current. Recall that CFT with c = 1 — 12«q < 
1 may be represented in terms of a free bosonic field f{z) with a back- 
ground charge 2a , see refs. [TUl HU E2] • The conformal weight of a state of 
coulomb charge a, or 2a — a, is h(a) = la(a — 2a ). The weight h r . s = 
[(rn — 4s) 2 — (k — 4) 2 ]/16/t corresponds to the charge a rjS = cto — \ct+ — |a_ 
with ct± the two screening charges. In the present case c = 1 — 6(n — A) 2 /4k 
and the correspondence is a_ = —2^2/k, a + = \J k/2 and 2a = a + + a_. 
In particular ai ;2 = — a_/2 and ao ; i/2 = «o i a_/4. The screening charges 
are such that the currents Q±(z) = exp(ia±(p(z)) have weight one. The op- 
erator $o in eq.(jBJ) is a linear combination of the primitive of the screening 
current Q- and the identity operator, ie: 

$ (z, z) = const'. 1 + 5te [const. / dw Q^(w) ]. 

Indeed this operator has conformal weight zero, satisfies the appropriate 
fusion rules and fulfills the charge conservation requirement which demands 
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that the sum of the coulomb charges of the operators involve in the correlation 
function minus the background charge belongs to the lattice generated by the 
screening charges. 



5 The limiting case k = 4. 

For k = 4 the SLE trace is a simple curve so that no point are swallowed and 
Pin = for all points. This case is marginal in the sense that the integral 
defining F(z) is only logarithmically divergent. By extension, we have: 



Pi(z, z) = -Sm 



TT 



log(tanh — ) 



121 



This satisfies the martingale equation (JJJ) for k = 4 and the appropriate 
boundary conditions: Pi(x G M+) = and Pi(x G R_) = 1. Contrary to the 
cases k < 4, it is discontinuous at the origin. On the upper boundary the 
distribution of the trace is given by: 

2 

Pi(iir + x) = 1 arctan(e :c/2 ), x G K. 

TT 

For k — 4, the Virasoro central charge is c = 1 and /ii ; 2 = 1/4 and 
ho-i/2 = 1/16. The probability fT2|) possesses a nice free field CFT inter- 
pretation. Central charge c = 1 corresponds to bosonic free field. Let us 
denote by X this field, hi-2 = 1/4 is the conformal weight of the boundary 
vertex operator Vi ;2 = cosX/\/2 which may be thought of as the boundary 
condition changing operator intertwining two boundary intervals on which 
two different Dirichlet boundary conditions are imposed, ho-i/2 — 1/16 is 
the dimension of the twist field a which is the boundary condition changing 
operator intertwining between Dirichlet and Neumann boundary conditions. 
Thus the probability Pi(z, z) is proportional to 

(X(z,z)) D;D;N = (a(x + )a(x-.)X(z,z)Vi ;2 (x )), 

where 'D;D;N' refers to Dirichlet boundary conditions on the lower boundary 
[x_,xo] and [xo,x + ], but with a discontinuity at xq (ie. two D-branes at 
finite distance), and Neumann boundary condition on the upper boundary 
(ie. a space filling brane). The fact Pi(z,z) satisfies the Dirichlet 
boundary conditions on the lower boundary is clear by construction but one 
may verify that it actually satisfies the Neumann boundary condition on the 
upper boundary. The fact that it is a harmonic function is then a consequence 
of the free field equation of X. 
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6 Boundary excursion probabilities (all ft's). 



The computations of previous sections yield the density probability of the 
hitting point on the upper boundary for k > 4. We now would like to show 
that this formula actually applies to any value of k. So we look for the 
probability P up (x), x G wr + 1R, that the SLE trace hits the upper boundary 
at a point u = m + u, u > x, on the right of x = m + x. By definition, this 
probability satisfies: 

P up (-oa) = 1, P up (+oo) = 0. 

To change gear we shall do the computation using conformal field theory 
techniques. Using the CFT martingales ©, or equivalently ©, we aim at 
proving that this probability is given by the correlation function 

(V , 0;l/2(^-)V , 0;l/2(x+)^o(^)V , l;2(^o)) ) 

where ipo(x) is a weight zero boundary conformal field inserted on the upper 
boundary. Recall that in the strip geometry x± = ±oo and Xq = 0. Due to 
CFT fusion rules, there are two possible choices for ip : either it is the iden- 
tity operator or it is the intertwiner between the Virasoro modules of weights 
hi-2 = and hvo = ^F^- The existence of two choices makes possible to 
fulfill the above two boundary conditions. These boundary condition code for 
the behavior of the functions as x approaches x_ or x + . In CFT correlation 
functions this behavior is governed by fusing ip with ipo-,1/2- in n on unitary 
CFT there is in general not much constraint on fusions. However, if the 
product of the operators ^o ; i/2^o is acting on the state created by ipi-2, 
then the null vector equation (— 2L_ 2 + ^L^^lu) = imposes constraints 
on the fusion. In the present case, taking into account that the out-state is 
created by tpip, we get: 

ipo-i/2(y)^o(x)\uj) c[V>o;i/2(2/) H + c (y - x) 2/K [i/j - 3 /2(y) H ]|w) 

where the dots refer to the descendant operators. The two fusion coefficients 
c and c depend on which operator ip we choose. The noticeable fact is that 
this operator product expansion is regular for all values of k : the coefficient 
in front of ipo-i/2 is constant and that in front of the V'o;3/2 vanishes as x — > y. 
This allows us to choose the operator ipo such that its fusion with -0 O; i/ 2 (a; + ) 
at point x + vanishes but that with -0o;i/2(^-) at point x_ tends to a constant. 
With this choice the correlation function satisfies: 

(V , 0;l/2(^-)^0;l/2(^+)V , o(^)V , l;2(^o)> 



J 1, if x —> x_; 
1 0, if x — > x + ; 
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By the martingale property (jSJ) or (jSJ), this correlation function but with x 
replaced by its image by f t , ie. by f t (x), is a martingale. It is such that, at 
large time, 

hm('0O;l/2(^-)'0O;l/2(^+)V'o(/t(^))V'l;2(^o)) = 1{ x on the left of 7sle}' 

t — >oo 

so that it projects on those events for which x is on the left of the SLE 
trace. Taking the expectation of this equation using the martingale property 
implies our claim that 

P up (x) = (V»o i i/2(aJ-)^V2(a?+)^o(^)V'i;2(jCo)>- (13) 

This CFT correlation function can be explicitly computed as the null 
vector relation (— 2L_ 2 + f-^-i)!^) — implies that it satisfies a second 
order differential equation. In the strip geometry the later reads: 

(^ + tanh(|)0 a )p tfP (z) = O, 

with x = m + x. Its solution is: 

For k > 4, it of course coincides with the formula we found in the previous 
sections, but the derivation we presented here is valid for any k. Surprisingly, 
the formula ()14|) shows no transition as a function of K, neither at k = 4 nor 
at k — 8. For k — 6, (1 — P up ) gives Cardy's crossing formula [H] for the 
existence of a cluster percolating from the boundary interval [x- , Xq] to the 
opposite boundary interval [x, x + ] in critical percolation. 

Why harmonic solutions to eq.(|7j) have a probabilistic relevance remains 
to be explained. It would have been legitimate to expect a simple formula 
also for the probability for a bulk point - and not only a boundary point - to 
be on the left of the SLE trace. Although this probability is again given by 
the bulk-to-boundary correlation function it is not harmonic for k < 4 
and we do not have a closed formula for it in this regime, as we do not have 
a clear understanding of this transition. Actually the same remark applies 
to the probability, meaningful for 4 < k < 8, that a bulk point be swallowed 
from the left by the SLE trace : it is not harmonic. 

7 Numerical simulations of Ising excursion 
probabilities. 

We have simulated the Ising ferromagnet on a square lattice strip of size 
L x 3L at the critical temperature (namely T c j J = 2/ log(l + V2), where J 
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is the ferromagnetic interaction). On the lower boundary spins, we applied 
a — J magnetic field on the left half of the boundary and a +J field on the 
right half of the boundary (this is equivalent to adding a raw of frozen spins 
below the lower raw of the strip). On the horizontal direction we applied 
antiperiodic boundary conditions, so that the minus phase on the left and 
the plus phase on the right would connect gracefully. This suppresses most 
finite size effects due to the finite length of the strip. The upper boundary 
was left free. 

To relax this system we used a cluster algorithm originally developed for 
spin glasses ^3] (using 64 configurations at one temperature) . This choice of 
algorithm was made because the code was available without much program- 
ming effort. More common techniques could have been used and would have 
probably lead to better running performance. 




Figure 2: A sample equilibrium configuration at L = 80, black and white 
squares respectively represent plus and minus spins. 

We have used different system sizes L = 10, 14, 20, 26, 40 and 80 and 
gathered for each size 320,000 independent samples (only 160,000 for L = 80). 
A sample configuration at L = 80 is shown on figure El For each sample, we 
measured the total lateral displacement of the interface. More precisely we 
followed the interface from the middle of the lower boundary up to the upper 
boundary and counted how many times it went to the left or to the right (the 
antiperiodic boundary may be crossed during this procedure). Note that the 
interface is somewhat ill defined on the square lattice: it can have branching 
points where it splits in two, in this case we chose one of the branch with 
probability 1/2. In some increasingly rare cases (less than 0.03% for the 
smaller size), it may also go around the system (through the antiperiodic 
boundary) and comes back to its starting point. In such cases we simply 
ignored the sample. 

We thus get the distribution of end point position x. The distribution 
Q(x r ) of x' = x/L is to be compared to the theoretical one computed previ- 
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Figure 3: Reduced distribution of the Ising interface end point (x' = x/L). 
The plain curve is the analytical result, data points are numerical results for 
L = 80. 
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Figure 4: Value of the maximum absolute difference between the numerical 
and analytical integrated distributions of the Ising interface end point for 
80 > L > 10. The straight line is just a guide to the eyes. 
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ously (pHj) (with n — 3) for a strip of width tt. In the simulation variables, it 
reads: 

Q( x >) = ^cosh(x'vr/2)~ 4/3 . (15) 

Where J is in this case 6 T(2/3) r(5/6)/v^F ~ 5.17422. The theoretical 
and numerical distributions (for L = 80) are shown on figure El They agree 
very well. To make a more precise statement, we considered the integrated 
distribution (for which the statistical error bars are much smaller) and com- 
puted the maximum 5 of the absolute difference between the theoretical and 
numerical integrated distributions. The results are shown on figure 0] This 
difference converges nicely to zero at large L. Moreover the finite size cor- 
rections seems to be of order 1/L. Note that the finite size corrections were 
expected to be at a least of order 1/L since the lattice spacing has a relative 
size of order 1/L. 



8 Conclusion. 

Besides chordal and radial SLEs, dipolar SLEs are the only possible simple 
SLEs on simply connected planar domains which satisfy a left-right reflection 
symmetry. As illustrated in previous sections, they provide a generalization 
of chordal SLEs with more structures since the hulls do not totally fill the 
domain for k > 4. The fact that the bulk visiting probabilities we computed 
are harmonic functions is clearly related to the discrete harmonic explorer. 
The latter was defined in ref.[13j to provide a discrete analogue of chordal 
SLE at k — 4. It clearly can be generalized to yield a discrete analogue 
of dipolar SLE at k = 4 by dealing with harmonic functions with 'D;D;N' 
boundary conditions to specify the probabilities of the excursion processes. 

The harmonicity property of the visiting probabilities for k > 4 was 
unexpected from a CFT point of view. The origin of the non-harmonicity 
transition - ie. the fact the probabilities for a point to be on the left (or on 
the right) of the curve stop to be harmonic for k < 4 - is not clear to us. We 
however noticed that for k > 4 the probabilities for a point to be swallowed 
from the left (or from the right) are also not harmonic functions. These two 
breakdowns of harmonicity seem to be related by duality k <-> 16/ k. 

Finally, numerical simulations of the Ising ferromagnet at criticality con- 
firms the analytical results presented here as well as the CFT mapping. 
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A Appendix: analytical details on the func- 
tion F(z). 

Here we gathered a few informations on the function 

du 

(sinhw/2) 4 / K ' 

We have to specify the analytical properties of (sinh z/2)~ i l K . It is such 
that arg(sinh z/2) G [0, 7r]. Hence, (sinh z/2)~ 4: ' K is real on the positive real 
axis, equals to e _i47r / K (sinh |x|/2) -4 / K on the negative real axis and equals to 
e _l27r//K (coshx/2)~ 4 / K on the upper boundary z = i-n + x. One may verify that 
it is equivalent to 2 4 / K e~ 2z l K around +00, while it behaves as 2 4//k e ( 22; - l47r )/ K 
around —00 . 

By analyticity, the contour of integration in the definition of F(z) can be 
choose arbitrary inside the strip §1 but starting at —00. As a consequence, 
comparing the integration along the real axis R and along ire + R we learn 
that F(0) = e~ H7r/K J and that 

F(+oo) = e~ i27r/K I = e~ i4n/K J + J, 

where / = dy (cosh y /2)~ 4 / K and J = J + °° <i?/(sinhy/2)~ 4 / K as in the 
text. This leads to / = 2Jcos(2tt/k) 

For z = a + ib, b G [0, tt], one may expand sinh z/2 as 

sinh z/2 = sinh a/2 cos b/2 + i cosh a/2 sin 6/2, 
so that (sinh z/2) = R~ i/K exp(-i46 z /K) with 

R z = I sinh z/2 1, 
tan# 2 = cotha/2 tan 6/2, 6> 2 G [0,7r]. 

The angle 6 a +ib decreases from (tt — b/2) to 6/2 as a varies from —00 to +00. 

The functions Shn[.F(z)] and Sm[e j2,r ^ F(z)] involved in eqs. (|lUlll|) . can 
be represented as: 

/X 
daR a+ib sin ( 4 #aW K )> 
-00 

/x 
daR^ sin((4^ +ib -2vr)// t ). 
-oo 

Other representations can be written using different integration contours. 
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